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Abstract 



We study incortimensurate fractional variational problems in terms of a generalized frac- 
tional integral with Lagrangians depending on classical derivatives and generalized fractional 
C^ I integrals and derivatives. We obtain necessary optimality conditions for the basic and isoperi- 

metric problems, transversality conditions for free boundary value problems, and a generalized 
Noether type theorem. 
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^ ■ 1 Introduction 

O 

Till recently, it was believed that Lagrangian and Hamiltonian mechanics were not valid in the 
presence of nonconservative forces such as friction [53]. In the last years, however, several ap- 
proaches have been investigated in order to find a Lagrangian or a Hamiltonian description for 
^ ' classes of dissipative (or dissipative- looking) systems t7i[8i|22j[271[37]. One possibility to have a 

^ . Lagrangian and a Hamiltonian formulation, for both conservative and nonconservative systems, 

was proposed by Fred Riewe in 1996 and consists in using fractional derivatives [351132] • Riewe's 
papers [381139] gave rise to a new and important research field, called the fractional calculus of 
variations [24| . Nowadays the subject is of strong interest, and many results of variational analysis 
were extended to the non- integer case (see, e.g., [2H51[^[T^I26 ,29,33]). Here we study problems of 
calculus of variations with generalized fractional operators ^1. .30..31j . Generalized fractional inte- 
grals are given as a linear combination of left and right fractional integrals with general kernels. 
Generalized fractional Riemann-Liouville and Caputo derivatives are defined as a composition 
of classical derivatives and generalized fractional integrals. In a first problem, we ask how to 
determine the extremizers of a functional defined by a generalized fractional integral involving n 
generalized fractional Caputo derivatives and n generalized fractional integrals. All these oper- 
ators have different (non-integer) orders. We obtain necessary optimality conditions, and in the 
case of free boundary values, also natural boundary conditions. Next, we derive Euler-Lagrange 
type equations for an extended isoperimetric problem and we obtain a Noether type theorem. 

The text is organized as follows. In Section[2]we give the definitions and main properties of the 
generalized fractional operators. We prove Euler-Lagrange equations for the fundamental gener- 
alized problem in Section [31 and natural boundary conditions for free boundary value problems in 
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Sectional Section[S]is devoted to the generalized isoperimetric problem and Section[B]to Noether's 
theorem. Finally, in Section [7] we present an application of our results to the damped harmonic 
oscillator. 

2 Preliminaries 

We start by defining the generalized fractional operators 1 . As particular cases, by choosing 
appropriate kernels, such operators are reduced to the standard fractional integrals and derivatives 
of fractional calculus (see, e.g., [2011211134) ). Throughout the text, a denotes a real number between 
zero and one. Following [Sj, we use round brackets for the arguments of functions, and square 
brackets for the arguments of operators. 

Definition 1 (The generalized fractional integral) . The operator Kp is given by 

X b 

K^ [f] (x) := K^ [t ^ fit)] {x)^p f k^{x, t)f{t)dt + q f k^{t, x)f{t)dt, 



where P — {a,x,b,p,q) is the parameter set (p- set for brevity), x e [a^b], p,q are real numbers, 
and ka{x,t) is a kernel which may depend on a. The operator Kp is referred as the operator K 
(K-op for simplicity) of order a and p-set P . 

Note that if we define 

pka{x, t) if i < X, 



^^^'^^' ^ qko.{t,x) iit>x, 
then the operator Kp can be written in the form 

K<^ [f] (x) =K^[t^ fit)] {x) - / G{x,t)f{t)dt. 

J a 

Thus, the generalized fractional integral is a Fredholm operator, one of the oldest and most 
respectable class of operators that arise in the theory of integral equations [16l|35] . 

Example 2. 1. Let ka{t - t) = r^{t - t)"-^ andO <a <1. If P = (a, t, 6, 1,0), then 

t 
1 /■<■ ^.-l 



a 

is the left Riemann-Liouville fractional integral of order a; if P — {a,t,b, 1,0), then 

b 

K^ifm = f^ /(^ - 0"- V(r)dr =: ,i?[/](0 
t 
is the right Riemann-Liouville fractional integral of order a. 

2. For ka{t -t) = r(a(Vr)) (^ " r)°(*-")-i and P ^ {a,t,b,l,0) 

t 

K^ifKt) = / ^^TTT^^^ - r)"^''^^-'f{r)dT =: aiT^^'-^i/KO 
J r(a(t,T) 

a 

is the left Riemann-Liouville fractional integral of variable order a{t,T), and for P 
(a, i, 6, 0,1) 

b 



t 
is the right Riemann-Liouville fractional integral of variable order a(t,T) \3^ . 



3. For < a < 1, ka{t^ r) = yCo) (^°S t) 7 '^^'^ ^ ^ ^'^^ ^' ^' -'^' ^)' ^^^ operator Kp reduces 
to the left Hadamard fractional integral \3^ . 

and for P = (a, t, b, 0, 1) operator Kp reduces to ttie right Hadamard fractional integral, 

4- Generalized fractional integrals can be also reduced to, e.g., Riesz, Katugampola or Kilbas 
fractional operators. Their definitions can be found in lU. 



Next results yield boundedness of the generalized fractional integral. 

Theorem 3 (cf. Example 6 of |16]). Let a £ (0,1) and P = {a,x,b,p,q). If ka is a square 
integrable function on the square A = [a,b] x [a,b], then Kp : L2{[a,b]) — >■ L2{[a,b]) is well 
defined, linear, and bounded operator. 

Theorem 4 (cf. [501151] ). Let ka G Li ([0, b — a]) be a difference kernel, that is, kdx, t) — ka{x—t). 
Then, Kp : Li {[a,b]) — >■ Li {[a,b]) is a well defined bounded and linear operator. 

Theorem 5 (cf. Theorem 2.4 of ^Slji)- Let P — {a,x,b,p,q). If ki^a is a difference kernel, 
/ci-a GLi([0,5-a]) and f €AC([a,b]), then Klr"[f] belongs to AC {[a, b]). 

The generalized fractional derivatives Ap and Bp are defined in terms of the generalized 
fractional integral K-op. 

Definition 6 (Generalized Riemann-Liouville fractional derivative) . Let P be a given parameter 
set and < a < 1. The operator Ap is defined by Ap := D o Kp~" , where D denotes the standard 
derivative operator, and is referred as the operator A (A-op) of order a and p-set P. 

Remark 7. Operator A is well-defined for all functions f such that Kp^^lf] is differentiate. 
Theorem assure us that the domain of A is nonempty. 

Definition 8 (Generalized Caputo fractional derivative). Let P be a given parameter set and 
a e (0,1). The operator Bp is defined by Bp := Kp~°' o D, where D denotes the standard 
derivative operator, and is referred as the operator B (B-op) of order a and p-set P. 

Remark 9. Operator B is well-defined for differentiate functions. 

Example 10. The standard Riemann-Liouville and Caputo fractional derivatives (see, e.g., \2(A 
\21[\84^ ) are easily obtained from the general kernel operators Ap and Bp, respectively. Let ka(t — 
^) = r(T^(^-^)~": «G (0>1)- IfP^{a,t,b,l,0), then 

t 

a 

is the standard left Riemann-Liouville fractional derivative of order a, while 

t 

'■ —a fi i^\j^ __. C 7-\Qr 



a 

is the standard left Caputo fractional derivative of order a; if P = {a, t, b, 0, 1), then 

b 

-A^Afm - -TXT^I /(- - ^ryir)dr =: ,Dn/](t) 

t 



is the standard right Riemann-Liouville fractional derivative of order a, while 

b 

^B-AfKt) = -Y^r^J(^'t)~"f'(^)dr =: ?DnfKt) 

t 

is the standard right Caputo fractional derivative of order a. 

The following theorems give integration by parts formulas for operators A, B and K. For 
detailed proofs we refer the reader to [301131] . 

Theorem 11. Let a £ (0,1), P = {a,t,b,p,q) , ka be a square-integrable function on A — 

[a, 6] X [a,b], and f,g € L2 {[a,b]). The generalized fractional integral Kf, satisfies the integration 

by parts formula 

b b 

I gix)K^ [f] ix)dx = I fix)K^, [g] {x)dx, (1) 

a a 

where P* =< a,t,b,q,p >. 

Theorem 12. Let a G (0,1), P = {a,t,b,p,q), and ka be a square integrable function on A — 
[a,b] X [a,b]. If functions f,Kp~"[g] € AC{[a,b]), then 



gix)Bf, [f] {x)dx = f{x)K'p-^ [g] {x)l - / f{x)A^p, [g] {x)dx, (2) 

J a 

where P* =< a,t,b,q,p >. 

Theorem 13. Let < a < 1, P — {a,x,b,p,q) , and ka be a difference kernel such that ka G 
Li[0, b — a]. If f € Li ([a, b]) and g € C ([a, 6]), then the operator Kp satisfies the integration by 
parts formula ([T]) . 

Theorem 14. Let a € (0, 1), P = (a, t, b,p, q) , and ka G Li ([0, b — a]) be a difference kernel. If 
functions f,g d AC{[a,b]), then formula ^ holds. 

For f = [/i, . . . , /jv] : [a, b] -^ M^, where TV e N, we put 



A"p [f] (x) 
B^ [f] (x) 
K^ [f] {x) 



[Af,[f,]{x),...,Af,[fN]ix)], 
[Bf.[f,]{x),...,B^[fN]{x)], 
[K^[f,]{x),...,K^[fN]{x)]. 



3 The generalized fundamental variational problem 

We consider the problem of finding a function y = [yi, . . . ,yN] that gives an extremum (minimum 
or maximum) to the functional 

J{y) ^K^[t^F {t, y(t), y'(i), B^J [y] (t), ...,B^pl [y] {t),Kl\ [y] (t), . . . , if^™ [y] {t))] (b) 

(3) 
subject to the boundary conditions 

y(a) = ya, y(&) = yb, (4) 

where a,Pi,jk S (0,1), P =< a, 6, 6, 1,0 >, Pi =< a,t,b,pi,qt >, and Rk =< a,t,b,rk,Sk >, 
i = 1, . . . ,n, k = 1, . . . ,m. For simplicity of notation, we introduce the operator {■}p p defined 
by 

{y}pn,H, (*) ■■- {t,y{t)yit),B^p^ [r ^ y(r)] {t),Kl [r ^ y(r)] (t)) , 



where 






(b^^...,b 









The operator Kp has kernel ka{x, t) and, for i = 1, . . . ,n and k = 1, . . . , m, operators Bp and 
if^'" have kernels hi-i3.{t,T) and h^^{t,T), respectively. In the sequel we assume that: 



(HI) the Lagrangian F e C^ {[a,h] x R^>^("+™+2);R); 



/3,7 



(H2) functions I? t h^ 9jv+jF{y}?,„%, Wfca(6, t) 



A' 



P^ 



T H> 



/3,7 



ka{b,T)d(,+i)M+jF {y}p2^p^ (r) 



^^- 



T H> fc„(6, r)5(„+i+fe)Ar+ji^ {y}pD,-Ri ('^) '^'^^ * '^ ^"(^' ^)^J-^ iylpJ^'.-R/ (^) ^^■'^ contin- 
uous on (a, 6), J = 2, . . . , A^ + 1, i = 1, . . . , n, fc = 1, . . . , m; 



1/9,7 



1/3,7 



(H3) functions 1 1-> Sat+jF {ylp^^^^^ {t)ka{b, t) and -?fp; t i-^ fca(&, T)d^^,+i^N+jF {yVp^^R^ (t) 
GAC([a,6]), j = 2,...,7V + l,^ = l,...,n; 

(H4) for i = 1, . . . , n, k = 1, . . . ,m, the kernels fcQ,(x, t), /ii-/3i (i, r) and h^,, {t, r) are such that we 
are able to use Theorems [Tl] [T^l [131 and/or HM 

Definition 15. A function y G C^ ([a, 6];R^) is said to be admissible for the fractional variational 
problem ^ -Q if functions Bp [y] and K'p [y], i = 1, . . . , n, fc = 1, . . . , to exist and are continuous 
on the interval [a, 6], and y satisfies the given boundary conditions Q. 

Theorem 16. If y is a solution to problem ©-(HI, then y satisfies the system of generalized 
Euler- Lagrange equations 

n 
k^{b,t)d,F {yf^2,Rr W -T.^% [^ ^ Ub,T)d^.+i)N+,F{yf^2,Rs W] W 



f2K]tf \r ^ k^{b,r)d^n+i+k)N+,F{yf^2M. (r)l (t) - - (On+.F {yf^^^Rr it)kc.{b,t)) = 



fc=i 



dt 



for allte{a,b),j = 2,...,N + l. 

Proof. The proof is analogous to that of 30 , Theorem 4.2]. 



(5) 



D 



4 Generalized free-boundary variational problem 

Assume now that in problem Q-O the boundary conditions ^ are substituted by 

y(a) is free and y(6) = yb- 



(6) 



Theorem 17. If y is a solution to the problem of extremizing functional (|3]) with ([6]) as the 
boundary conditions, then y satisfies the system of Euler-Lagrange equations ^ . Moreover, the 
extra system of natural boundary conditions 

n 



j = 2,...,7V+ 1, holds. 

Proof. The proof is analogous to that of JSOl, Theorem 5.1]. 



D 



5 Generalized isoperimetric problem 

Let ^ e M. Among all functions y : [a, b] — >■ E^ satisfying the boundary conditions 

y(a) = ya, y{b)=yb, 

and an isoperimetric constraint of the form 

we look for those that extremize (i.e., minimize or maximize) the functional 

J{y)=K^\F{yfplJib). 



(8) 



(9) 



(10) 



For i — 1, . . . ,n, k = 1, . . . ,m operators Kf,, Bp, and K'^ , as well as function F, are the same as 
in problem ©-(Hj). Moreover, we assume that functional (HJ satisfies hypotheses (H1)-(H4). 

Definition 18. A function y : [a,b] — s> R^ is said to be admissible for problem ((8|)- (|10l) if 
functions Bp, [y] and K]^ [y], i = 1, . . . ,n, k = 1, . . . ,m, exist and are continuous on [a, b], and y 
satisfies the given boundary conditions ^ and the isoperimetric constraint ^. 

Definition 19. An admissible function y ^ C^ ([a, &],R^) is said to be an extremal for I if it 
satisfies the system of Euler-Lagrange equations ([5]) associated with functional in . 



Theorem 20. If y is a solution to the isoperimetric problem 
X, then there exists a real constant A such that 



TU)) and is not an extremal for 



k^{b,t)d,H{yf^2,Rr (t)+J2Kli k ^ Ub,r)d(^,,+,+k)N+jH{yyp2^p^ (r)j (t) 



k=l 



-J2A^, \t ^ k^ib,T)d^,+r)N+,H{yf^lp^ (r)! (i) - - (o.+nH {yfp^p^ it)k^{b,t)) = 



for all t G (a, 6), j = 2, . . . ,N + I, where H{t, y, u, v, w) — F{t, y, u, v, w) — AG(i, y, u, v, w), 
P* =< a,t,b,qi,pi >, and Rl =< a,t,b,Sk,rk >■ 



Proof. The proof is analogous to that of ^30j Theorem 6.3]. 



□ 



6 Generalized fractional Noether's theorem 

Emmy Noether's theorem on extremal functionals, establishing that certain symmetries imply 
conservation laws (constants of motion) , has been called "the most important theorem in physics 
since the Pythagorean theorem" . For a recent account of Noether's theorem and possible appli- 
cations in physics, from many different points of view, we refer the reader to [5S]. Formulations 
in the more general context of optimal control can be found in [13130] . Conservation laws appear 
naturally in closed systems. In presence of non-conservative or dissipative forces, the constants 
of motion are broken and Noether's classical theorem ceases to be valid. It is still possible, how- 
ever, to obtain Noether type theorems that cover both conservative and non-conservative cases. 
Roughly speaking, one can prove that Noether's conservation laws are still valid if a new term, 
involving the non-conservative forces, is added to the standard constants of motion [TU]. The first 
Noether theorem for the fractional calculus of variations was obtained in 2007 [TT]. Since then, 
the subject attracted a lot of attention. The state of the art is given in the book [24]. Here we 
obtain a Noether's theorem for generalized fractional variational problems. 



Definition 21. We say that the functional ^ is invariant under an e-parameter group of in- 
finitesimal transformations 

y(i)=y(i)+e^(t,y(t))+o(£) (11) 

if for any subinterval [ta,tb\ ^ [a,^] one has 

Kt [t ^F{yfp2,B, (O] (4) = ^1 [t ^F{yfp2.n, (O] {hi 
where P = {ta,tb,tb, 1,0). 



(12) 



Tiieorem 22. If functional ^ is invariant under an e-parameter group of infinitesimal trans- 
formations, then 



N+l 



E r^-^Mp.%. (0 •O-i(i.yW) + a^+,^{y}^;^%^^ (t) . e,_i(t,y(t)) 

j=2 \ 

n 

+E%+i)^+.^{y}p;.^,«. w • s^:[t ^ e,-i(r,y(r))](i) 

m \ 

+ E d^^+,+,-,M+,F{yfp2Ms W • if;i:[T ^ e,-i(r,y(r))](t) ^ 0. (13) 



fe=i 



Proof. Since, by hypothesis, condition ([T2l) is satisfied for any subinterval [ta,tb\ C [a,b], we have 

F{yfpluAt)-F{yfplHAt)- (14) 

Differentiating (1141) with respect to e, then putting e = 0, and applying definitions and properties 
of generalized fractional operators, we obtain p^ . D 



In order to state the Noether theorem in a compact form, we introduce the following operators: 

1 



Dp[/,ff]W 
Ip[/,5](t): 



ka{b,t 

-1 



■fit)-A'^p.[g]it)+git)-B^[f]it), 



fit)-K^4g]it)+g{t)-K^[f]it), 



ka\b,t) 
where P* denotes the dual p-set of P, that is, if P = (a, i, h,p, q), then P* = (a, t, 6, q,p). 



(15) 
(16) 



Tiieorem 23 (Generalized fractional Noether's theorem). If functional ([3]) is invariant under an 
e-parameter group of infinitesimal transformations (jlip . then 



N+l / n 



E E^^p: r "^ ^,-i(r,y(T)),T K^ A;„(6,r)%+i)^+,P{y}^;^%^ (t)J (f) 

+ E^Pfc f^ ^ 6-i(^,y(T)),T ^ fca(6,r)9(„+i+fe)Ar+jP{y}^;^%^ (r)l (t) 



fe=i 



1 d 



+ ^ (O-i (t,y(t)) • dN+,F{yfp2,B., {t))H,-iitMt)ydN+,F{yfp2,,^ (^ty——-k^{b,t) ^ 



(17) 



for any generalized fractional extremal y of J and for all t G (a, 6). 



Proof. By Theorem 1161 we have 

n 



i^l 



J2k]^ [t ^ k^{b,r)d^n+i+k)N+,F{yf^lj,^ (r)] (t) + - (9^^+,^ {y}^;^%^ (0^.(6,0 



fe=i 



(18) 



for all i e (a, 6), j = 2, . . . , Af + 1. Substituting ([TH) into ([12]), we obtain 

JV+l 



E 

i=2 



1 



ka{h,t) 



\i=l 



^j-i{tMt)) E4'* ^ ^ fc„(&,r)%+i)^+,i^{y}^;J,^^ (r) (i) 



Y,K]{1 [t ^ k^{b,T)d^,,+,+,^^+,F{yyj;l^^ (r)] (i) + - (a^+,F {y}^;^^^^^ (i)fca(fe,i)) 



k=l 



dN+,F {yfplj,^ it) ■ -e,-i {t, yit)) + ^ 9(,+i)A.+,F {y}^;^^^^^ (i) • B^; [r .^ O-i (^, y(^))] (*) 



J29in+^+m+,F{y}po.Rr W • A'S;[r ^ e,-i(r,y(r))](i) 



fc=i 



Finally, we arrive to (HH) by (HSJ and ([T6l) . 

Example 24. Lei P = (a, i, 6,p, g). Consider the following problem: 



D 



(19) 



(20) 



J[y]^ / P(i,B^[y](t))dt^min 

J a 

y{a) = ya , y{b) = 2/6, 

ararf transformations 

y{t) = y(t) + £c + o(e), 

where c is a constant. For any [ta,ti,] C [a, 6] we have 

[' F{t,B^[y]{t))dt= [" F{t,B^mt))dt. 

Jta Jta 

Therefore, J^[y] is invariant under (|20p and Theorem \23\ asserts that 

AUr^d,F{T,Bf,[y]{Tmt)^0 (21) 

along any generalized fractional extremal y. Notice that equation (I21|l can be written in the form 

^^{K^,[r^d,FiT,B'^Ay]irmt)) = 0. 

In analogy with the classical approach, quantity Kpt[T — > d2F [t, Bp[y](T))]{t) is called a gener- 
alized fractional constant of motion. 



7 Applications to Physics 

If the functional ([3]) does not depend on B-ops and iiT-ops, then Theorem [TBI gives the following 
result: if y is a solution to the problem of extremizing 



J{y)= / L{t,y{t),y'{t))k^{b,t)dt 



(22) 



subject to y(a) — ya and y(6) = y^, where a G (0, 1), then 

d,L{t,y{t),y'{t)) j^dN+,L{t,y{t),y'{t)) = ^-1^ ■ ^^Ub,t)dN+jL {t,y{t),y' (t)) , (23) 



j ~ 2, . . . ,N + 1. In addition, if we assume that functional (j22p is invariant under transformations 
dTTI) . then Noether's theorem yields that 



(^(e,-i(i,yW)-9jv+,i(i,y(i),y'(t)))+e,-i(t,y(t))-aAr+,i(t,y(t),y'(i)).-^-^ = 0, 



E 

along any extremal of ([2^ . Let us consider kernel ka{h,t) — e"^^~*) and the Lagrangian for a 
three dimensional system: 

L (y, y) = -m (yi^ + ya^ + ^32) _ v{y), 

where V{y) is the potential energy and m stands for the mass. Observe that an explicitly time 
dependent integrand L = c°'^''^*''>L of functional ((22|) is known in the literature as the Bateman- 
Caldirola-Kanai (BCK) Lagrangian of a quantum dissipative system [TSlIll]- But in our case the 
Lagrangian of the system is L and not e"^''~*^L. The Euler-Lagrange equations (f23|) give the 
following system of second order ordinary differential equations: 

'y,{t)-am{t) = -7^d^V{y{t)) 

< y2it)-ay2{t)^-^d2V{y{t)) 

Mt)-^y3(t)^-T^93V{y{t)). 

If 7 := —a, then 

1 dV 

j/»+7y» + -T^ = 0, (24) 

i — 1,2,3, which are equations for the damped motion of a three-dimensional particle under the 



action of a force 



dv dv av 
into equations for a damped simple harmonic oscillator 



see, e.g., [I7|). Choosing V := k ^^ ^^ ^^ , we can transform 

iSC 

1 = 1,2, 3, with aj2 — A Now, let us consider the following Lagrangian: 

L (y, y) = 2™ (y'l^ + y'2^ + y^^) ~ ^gyl- (25) 

We see at once that the Lagrangian (1251) is invariant under the transformation 

yi^yi + £, m = 1/2, ys = 2/3- 

In this case Noether's theorem gives 

— (myi) ^amyi. (26) 

at 

If a = 0, then there is no friction and (j26p yields the classical conservation of linear momentum 

Pi — myi = const. Observe that the generalized momentum conjugate to yi is Pi — ^ — myi, 

i = 1,2, 3. This is not the case for the the BCK Lagrangian [T3], where the canonical momentum 

for j/j is Pi = e°''^'^~*^myi, i = 1,2,3, that is different from the kinetic momentum. Now, let us 

suppose that L is variationally invariant under the transformation 

yi = yi cos e + 2/2 sin e, -1)2 = -yi sin e + y2 cos e, ys ^ ya- 

Then ^1 = j/2, ^2 — —yi and ^3 = 0. For this case Noether's theorem yields 

-^{myiy2 - myiy2) - am{yiy2 ~ ym) = 0. (27) 



Note that for a — relation (j27p gives the standard conservation law piy2 — P2yi — const yielded 
by the classical Noether's theorem [HI Section 9.3]. 
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